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1. Abstract 

In this work, the performance of Tunnel Field-Effect 
Transistors (TFETs) based on two-dimensional Transi¬ 
tion Metal Dichalcogenide (TMD) materials is investi¬ 
gated by atomistic quantum transport simulations. One 
of the major challenges of TFETs is their low ON- 
currents. 2D material based TFETs can have tight gate 
control and high electric fields at the tunnel junction, 
and can in principle generate high ON-currents along 
with a sub-threshold swing smaller than 60 mV/dec. 
Our simulations reveal that high performance TMD 
TFETs, not only require good gate control, but also 
rely on the choice of the right channel material with 
optimum band gap, effective mass and source/drain dop¬ 
ing level. Unlike previous works, a full band atomistic 
tight binding method is used self-consistently with 3D 
Poisson equation to simulate ballistic quantum transport 
in these devices. The effect of the choice of TMD 
material on the performance of the device and its transfer 
characteristics are discussed. Moreover, the criteria for 
high ON-currents are explained with a simple analytic 
model, showing the related fundamental factors. Finally, 
the subthreshold swing and energy-delay of these TFETs 
are compared with conventional CMOS devices. 

II. Introduction 

Power consumption is one of the main challenges 
for future electronics. In this regard, tunnel FETs are 
among the most promising candidates for future in¬ 
tegrated circuits (ICs) due to their small subthreshold 
swing (SS) and low OFF-current HD, O- Having small 
SS and OFF-current reduces both static and dynamic 
power consumption of the ICs O. 


One of the major drawbacks of tunnel FETs is their 
low ON-current. The current in TFETs is the result 
of band-to-band tunneling (BTBT) of the carriers. The 
tunneling probability is usually much smaller than 1 and 
as a result the ON-currents of TFETs are much smaller 
than those of conventional FETs. However, the tunneling 
probability can increase significantly if the electric field 
at the tunneling region is high enough. Atomically thin 
2D devices are very interesting in this context for TFET 
applications, analogous to nanotubes HI, O, due to 
the tight gate control over the channel which results 
in high electric fields at the tunnel junction. Since the 
BTBT transmission probability depends on the electric 
field exponentially, high ON-currents are expected in 
2D TFETs. However, tight gate control is not the only 
player in determining ON-currents and there are other 
critical factors, such as band gap, effective mass (m*), 
and doping concentration. 

The well known scaling length theory O can be used 
to quantify the effect of gate control on the electric 
field at the tunnel junction ITl. This theory provides a 
simple analytic way to understand how various device 
parameters affect the spatial variation of the potential 
along the channel ||6l described by a modified ID 
Poisson equation, 
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where V and A are the electrostatic potential and the 
natural scaling/decay length of the potential, respectively. 
In double gated FETs, A is given by lO 
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where Cch, and Cqx are the dielectric constants of channel 
and oxide respectively, while tch, and tox are their 






thicknesses. According to equation reducing the 
channel thickness reduces the natural scaling length of 
the potential which results in higher electric fields and 
ON-currents. 

Besides advantageous 2D electrostatics, there are other 
incentives for using 2D materials for future electronics 
ii. For example, thinning the 3D material based FETs 
(which is required in transistor scaling) increases the 
effect of surface roughness on carrier transport which 
leads to lower mobility nni. However, 2D materials have 
relatively weak inter-layer bonds and can be exfoliated 
easily without surface roughness lITTIl . (121, CSI, (FTl - 
Moreover, there are no dangling bonds in 2D TMDs 
unlike thinned 3D materials ca, ca, ca. Another 
advantage of 2D materials is that thinning the material 
does not increase the band gap of the material as much 
as it does in 3D materials ca, ca. This is again 
due to the weaker coupling between stacked layers in 
2D materials. For example, thinning In As nanowires to 
achieve better gate control can increase the band gap 
more than 100% d, while the increase in band gap 
from bulk to monolayer is much smaller in 2D materials 
f\M . This is useful in the TFETs as a larger badgap 
results in a lower ON-current. Moreover, mono-layer 2D 
materials usually have small dielectric constants ||20l, 
which can also increase the ON-current and reduce the 
drain induced barrier lowering (DIBL) in TFETs. 

Previously, TMD TFETs have been simulated without 
solving the Poisson and the transport equations self- 
consistently (H, O. In these works, the electric field 
at the source-channel junction has been approximated 
assuming that the ID Poisson equation (equation Q) 
is accurate throughout the entire source-to-channel junc¬ 
tion. Although A is a critical factor in determining the 
electric field, it is not the only factor. Other factors such 
as the depletion width in the source region are important 
too. The non self-consistent results predict very large 
ON-currents for all TMD TEETs irrespective of the 
channel material. In this paper, the correctness of these 
assumptions for TMD TEETs has been investigated by 
solving the 3D Poisson equation self-consistently with 
full band quantum transport, and the factors limiting the 
ON-currents are elucidated. 

The band structure and electronic properties such as 
band gap, m*, and dielectric constant of TMD materials 
depend on the number of layers. Consequently, devices 
with different number of layers show different charac¬ 
teristics. Although increasing the thickness reduces the 
band gap, it increases A and decreases the electric field 
at the source-to-channel junction. Moreover, multi-layer 
TMDs usually exhibit indirect band gaps which implies 


that phonons need to be involved in the BTBT process. 
As a result, it is favorable to use mono-layer TMDs 
and instead explore different materials to obtain small 
values of A, direct band gap, and m* simultaneously. 
Since TMDs form a general class of materials of the 
form MX 2 , where M is a transition metal (Mo, W, etc.) 
and X is a Chalcogenide (Te, Se, S), a variety of material 
parameters can be accessed by the correct choice of 
material. The field of 2D materials is still at its infancy 
as novel materials are being discovered (23l, (24l, which 
opens up opportunities for TEET designs. Accordingly, 
in this paper, only mono-layers of a set of more common 
TMDs, whose critical parameters span the design space 
of TEETs, are studied. 

III. Simulation method 

The TMD Hamiltonian is represented by an sp^d^ 2nd 
nearest neighbor tight-binding (TB) model with spin- 
orbit interaction. The Slater-Koster TB parameters 
are optimized based on first principles bandstructures 
obtained from density functional theory (DET) with the 
generalized gradient approximation (GGA). DET-GGA 
has been shown to provide band gaps and effective 
masses in TMDs comparable to experimental measure¬ 
ments (261 . The motivation for using a DET guided TB 
model is that a realistically extended device size can be 
simulated at ease compared to computationally expensive 
and size limited ab-initio methods Bandgaps and 
effective masses of mono-layer M 0 S 2 , WSe 2 , MoTe 2 , 
and WTe 2 obtained from our TB model are listed in 
table |T| The TB parameters are general and capture the 
bandstructure of both bulk and monolayer TMDs. As an 
example, the TB bandstructure of monolayer WTe 2 is 
shown in Eig. 

In this work, self-consistent Poisson-QTBM (Quan¬ 
tum Transmitting Boundary Method ||28l) methodology 
has been used within the tight binding description. The 
QTBM method is equivalent to the well known non¬ 
equilibrium Green’s function (NEGE) approach without 
scattering, but is a more computationally efficient imple¬ 
mentation II 29 II . In this method, the Schroedinger equa¬ 
tion with open boundaries is solved using the following 
equation 

{El-H- E)^s/d = Ss/D (3) 

where E, /, H, and X are energy, identity matrix, 
device Hamiltonian, and total self energy due to open 
boundaries and 7 /^ and S are the wave function and 
a career injection term from either source (subscript 
S) or drain (subscript D) (23. The electron and hole 
carrier density and current can be obtained from the 
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wave function 7 /^ 1^ . Since in-plane and out-of-plane 
dielectric constants and of TMD materials 

are different, the Poisson equation reads as follows if 
considering the z direction to be along the out-of-plane 
direction (or c-axis) of the TMDs 
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where V and p are the electrostatic potential and total 
charge, respectively. The dielectric constant values 
and for TMD materials are taken from ab-initio 
studies 1 ^ . 1211 and are listed in table |T| In this 
work, transport simulations have been performed with 
the nanodevice simulation tool NEMOS m, ED. 


of 0.5V is used unless mentioned otherwise. Equivalent 
oxide thickness (EOT) of top and bottom oxides is set to 
0.43 nm to be consistent with International Technology 
Roadmap for Semiconductors (ITRS) projections for 
2027 03. 


Gate 


p 


Top oxide 

N 

J 




Source 


Bot. oxide 

:]- 

Drain 


Vds 


TABLE I: Band gap (Eg), electron and hole effective 
masses (m* and and in-plane and out-of-plane rela¬ 
tive dielectric constants (e*^ and e^^^) of TMD materials. 
All the parameters are obtained from TB bandstructure 
of TMD monolayers with the exception of the dielectric 
constants which are taken from ab-initio studies 1201 . 

ED. 


Parameters 

M0S2 

WSe2 

MoTe2 

WTe2 

Eg [eV] 

1.68 

1.56 

1.085 

0.75 

ml [mo] 

0.52 

0.36 

0.57 

0.37 

ml [mol 

0.64 

0.5 

0.75 

0.3 

An 

4.2 

4.5 

8 

5.7 

f-OUt 

2.8 

2.9 

4.4 

3.3 



Eig. 1: Bandstructure of a mono-layer WTe 2 . 


IV. Results and discussion 

All simulated TMD TEET devices assume a structure 
as shown in Eig. and have channel and source/drain 
lengths of 15nm and lOnm, respectively. A doping level 
of le20 cm“^ is assumed in the source and drain regions 
which seems feasible by molecular doping of source and 
drain contact regions ISD . A source-drain voltage Yds 


Eig. 2: Physical structure of a mono-layer TMD TEET 
with channel and source/drain lengths of 15nm and 
lOnm, respectively. Source and drain regions have a 
doping level of le20 cm“^ and EOT is fixed to 0.43 
nm. 


A. Transfer characteristics of TMD TFETs 

Eig. shows the transfer characteristics of the TMD 
TEETs with OEE-current fixed to 1 nA/pm at 0 gate 
voltage. It is worthwhile to notice that in TEETs lower 
OEE-currents can be readily achieved without losing 
too much ON-current. This is because the slope of the 
Ids — Vcate curve is very steep in the low current 
regime. As a result, a much lower OEE-current can 
be obtained with a very small reduction of Vcate- For 
example, an OEE-current of 0.1 nA/pm (10 times 
smaller than before) can be obtained with the ON-current 
just being approximately 5-10% smaller. 

The simulation results show that WTe 2 TEETs can 
provide highest performance in terms of ON-current and 
SS in comparison to the other TEETs. Since WTe 2 has 
the smallest band gap and effective mass compared to the 
other TMDs, its ON-current is significantly higher. No¬ 
tice that despite the fact that MoTe 2 has a smaller band 
gap than WSe 2 , it shows a smaller current. The values of 
ON-current (the current at VGs=yDS=yDD=^-^Y), band 
gap and reduced effective masses (1/m* = 1/m* -1- 
l/ml) of these TMDs are listed in table Although 
M 0 S 2 has a high effective mass for tunneling applica¬ 
tions, it is ideal for ultra-scaled MOSEET applications 
where a high effective mass can suppress the source 
to drain tunneling |[34l . ESI. To understand the origin 
of the difference between ON-currents of these devices, 
one needs to consider several factors: electron and hole 
effective masses, band gap, body thickness, EOT of 
oxide, source-to-channel potential difference, source and 
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mining the current in TFETs is 




Fig. 3: Transfer characteristics of TMD TFETs with 

IqFF — ^TiA/(ITY l. 

TABLE II: ON-current Iqn [/iA//im], band gap Eg 
[eV], reduced effective masses m* [mo], natural scaling 
length A [nm], effective natural scaling length A [nm] 
from simulation and r] values [unitless] of TMD TFETs 


Material 

Ion 

Eg 

m* 

A 

A 

1 

WTe2 

127 

0.75 

0.17 

0.45 

2.45 

3.15 

WSe2 

4.6 

1.56 

0.21 

0.41 

2.5 

5.2 

MoTe2 

2.3 

1.08 

0.32 

0.5 

2.7 

5.8 

M 0 S 2 

0.3 

1.68 

0.29 

0.38 

2.5 

6.3 


drain doping level. These factors are considered in the 
analytic equation for the current in the ON-state of a 
TFET (361: 


I oc exp 


^ 3qhE J 


(5) 


where q, Ti, and m* are the charge of an electron, the 
reduced Plank constant, and the reduced effective mass, 
respectively. E is the electric field at the source-channel 
junction and can be approximated by (0^ — 
where A is the effective natural scaling length of the 
potential at the source-channel junction and q{(t)s — (t>ch) 
is the source-channel potential difference. At threshold 
voltage, q{(j)s — (j^ch) equals Eg. For small overdrive 
voltages (Vcs — yTh)^ the current can be simplified using 
q{<l>S - (t>Ch) ~ Eg to 


/ (X exp 


-4A^2m*Eg\ 

sn J 
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This equation shows the fundamental factor in deter- 


Ay^m^Eg 
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The factor r] depends on the device design (EOT, body 
thickness, and doping) through A and band structure 
of channel through m* and Eg. Note that for small 
overdrive voltages, the band gap and the effective mass 
have an equal impact on the current. For example, 
despite the fact that MoTe 2 has a smaller band gap 
compared with WSe 2 , it has a larger reduced effective 
mass which can partly compensate for the reduction in 
the band gap. 

The parameter A is composed of two components: 
the depletion width in the source region (Wd) and 
the natural scaling length (A) in the gated region. The 
difference between this analysis and previous works (H, 
0 , (33 is that the effective natural scaling length A is 
calculated from a self-consistent potential profile which 
includes the effect of both Wd and A accurately. Ideally, 
the dielectric constant of the channel should be small 
to make A as small as possible. Wd and A depend 
on the in-plane and out-of-plane dielectric constants, 
respectively. The difference between the values of A 
and A, listed in table [n| shows the importance of Wd- 
According to DFT simulations, MoTe 2 has the higher 
dielectric constant if compared to WSe 2 (^ which 
further lowers its ON-current. The factor rj for these 
materials is shown in table |I^ The larger the r], the 
smaller the current. By comparing the r] values, it can 
be seen that WSe 2 provides higher currents close to the 
threshold voltage compared to MoTe 2 TFETs in spite of 
its higher bandgap. Moreover, it is expected that WTe 2 
TFET produces the highest ON-current. Note that by 
employing equation ([^, the ON-current ratios between 
TMD FETs from different materials can be reproduced 
within reasonable accuracy when using the rj values from 
table ini 


Erom Eig. one can conclude that WTe 2 TEETs 
show promising ON-currents and a steep SS compared 
to other TMD TEETs. It is possible to further increase 
the ON-current of WTe 2 TEET by increasing the doping 
concentration at the source and drain regions. Increasing 
the doping decreases the depletion width of the source- 
to-channel interface, which also decreases r] according 
to equation (|7]). Consequently, the ON-current increases. 
Eig. 1^ shows a comparison between the ON-currents of 
WTe 2 TEETs with source/drain doping levels of le20 




cm“" and 2e20 cm”"^. In the case of 2e20 cm”"^ doping 
level, the ON-current increases to 350 [lAj[im which is 
a very high ON-current if compared with other TEETs. 
It is important to notice that the electric field at the 
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source-to-channel interface depends on both the natural 
scaling length A and the depletion width Wd- a. result, 
reduction of A is not sufficient for high ON-currents. A 
smaller Wd is also needed. In this regard, the doping of 
the source and drain region should be designed carefully, 
whether this doping is chemically or electrically induced. 
Fig. shows the band profiles for these two different 
doping concentrations. It is apparent that the higher 
doping concentration shows higher electric field due to 
the smaller Wd- Notice that the drain region is usually 
doped less than the source region to suppress the p- 
branch in the I-V characteristics and improve the OFF- 
state performance of the TFET |[39ll . However, due to a 
relatively large band gap and effective mass of TMD 
materials, the minimum achievable current in TMD 
TFETs is always very small and below 1 nAj[im even 
for a drain doping level of 2e20 cm“^. Consequently, 
lowering the drain doping does not have any significant 
impact on the transfer characteristics of TMD TFETs. 



Eig. 4: a) Transfer characteristics and b) band diagrams 
of WTe 2 with doping levels of le20 cm“^ and 2e20 
cm“^. 


B. C-V and DIBL 

Eig. shows the total gate capacitance versus the 
gate voltage (C-V) for the TMD TEETs under investi¬ 
gation. There are two major factors which determine the 
C-V characteristics: the quantum capacitance, and the 
threshold voltage. The total gate capacitance (Cq) can 
be modeled as a series of oxide capacitance (Cqx), and 
quantum capacitance (Cq). The quantum capacitance 
is proportional to the density of states and m* of the 
material 1401 . ED. Materials with larger m* have larger 
Cq which translates into larger Cq. On the other hand, 
the threshold voltage determines the voltage where the 
charge starts to appear in the channel. Since the off- 
currents of TEETs have been fixed here, the threshold 


voltage of one material is different from the other. 
Clearly, WTe 2 is showing the lowest Cq values, another 
benefit when it comes to benchmarking of various TMD 
devices as discussed below. 



Eig. 5: C-Wcate of WTe 2 , WSe 2 , and MoTe 2 TEETs. 

Drain induced barrier lowering (DIBL) is one of 
the most important short-channel effects for ultra-scaled 
transistors. DIBL causes a reduction of the threshold 
voltage by the drain voltage, and it is one of the com¬ 
monly used criteria for indicating short-channel behavior 
1^ . In TEETs, the active region in which tunneling oc¬ 
curs is right at the source-to-channel interface, far from 
the drain contact. Moreover, in 2D TEETs, gate control 
is stronger due to a thin channel which suppresses short 
channel effects. Consequently, short channel effects are 
less significant in comparison with conventional EETs 
(n-i-n or p-i-p doped transistors). The numerical values 
of DIBL for TMD TEETs obtained in this work are listed 
in table III. Note that DIBL values of WSe 2 and WTe 2 
TEETs are substantially smaller than the reported 80 
mVA^ DIBL of ultra-scaled MOSEETs El, El. DIBL 
is calculated at the current level of InA/iim from the 
change in the threshold voltage caused by varying Yds 
from O.IV to 0.5V. 

]jj^l ^ ^ThiVps = 0-5) - VThiVps = 0-1) .g. 

0.5-0.1 

Among these TMD TEETs, WSe 2 and WTe 2 show 
a smaller DIBL compared to MoTe 2 . This is because 
WSe 2 and WTe 2 have a smaller in-plane dielectric 
constant compared to MoTe 2 which reduces the electric 
field penetration from drain, and hence suppresses short 
channel effects. 
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TABLE III: DIBL values for TMD TFETs 


Material 

WTe2 

WSe2 

MoTe2 

DIBL [mVA^] 

25 

20 

67 


C. Subthreshold swing and energy-delay product 

The main idea behind a TFET is to achieve a steep 
sub threshold swing below the Boltzmann limit of 60 
mV/dec. It is important to compare the steep devices 
in a generic way. For example, the average subthreshold 
swing does not provide information about the current 
range in which the I-V is steep. Recently, a method for 
benchmarking steep devices has been proposed which 
gives insight into the local steepness of the I-V curve 
ll43l . In this method, the subthreshold swing is plotted 
against the drain-to-source current density of the device. 


Fig. 6a shows how mono-layer TMD TFETs fit in this 
picture. Notice that the current on the horizontal axis is 
not the ON-current but the current where the SS has been 
calculated. The devices with the best performance will 
ideally be placed on the lower right comer of the figure, 
where the SS is small and the drain-to-source current is 
large. 

Fig. shows again that mono-layer WTe 2 is the 
most promising candidate for TFET applications among 
the TMD materials considered in this study. Fig. [6b| 
shows the energy versus the intrinsic delay 

(pGateyDDlIoN) of a WTe 2 TFET versus ultra-scaled 
silicon MOSFETs. The WTe 2 TFET shows a smaller 
intrinsic energy-delay product value when compared to 
MOSFETs with the same channel length. In terms of the 
anticipated circuit performance of a 32-bit adder based 
on TMD TFETs, the performance specs discussed in this 
article for WTe 2 devices for Vdi)= 0.5V result in the sim¬ 
ilar energy-delay product as reported by Nikonov et al. 
1441 for VDi:)=0.25V. The actual energy and delay values 
were calculated to be 13fJ and 2450ps respectively using 
the same code as in (44). 


V. Conclusion 
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Fig. 6: a) SS-Ids for different TMD TFETs b) Energy- 
delay of WTe 2 TFET with a doping level of 2e20 cm-3 
versus Si MOSFETs from ITRS 1331 . |4^ . The dashed 
lines depict constant energy-delay products. 


effects: DIBL and SS values are significantly lower than 
for Si MOSFETs (by a factor of about 1/3). Moreover, 
the energy-delay product of the optimized WTe 2 TFET 
is lower than that of an ultra-scaled Si MOSFET. Our 
simulations show that 2D materials with lower band gaps 
(0.5-0.7eV) and effective masses are more suitable for 
high performance TFETs. 
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In this work, the performance of various Transition 
Metal Dichalcogenide materials (M 0 S 2 , WSe 2 , MoTe 2 , 
and WTe 2 ) TFETs has been investigated through self- 
consistent atomistic simulations. It has been shown that 
an atomically thin channel alone is not sufficient for 
high performance TFETs; especially to achieve high 
ON-currents, the choice of channel material and device 
design are critical as well. According to our analy¬ 
sis, WTe 2 is the most promising TMD in this study 
for TFET applications with high ON-currents of 350 
pA/pm. TMD TFETs exhibit reduced short channel 
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